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Abstract 

The energy spectrum of Mkn501 in the TeV energy regime, as measured by the 
HEGRA (High Energy Gamma Ray Astronomy) Gerenkov telescopes during its 
low state in 1995/96 and during a fraction of the 1997 outburst in the TeV energy 
regime, is shown to place stringent upper limits on the still unknown infrared photon 
density in the energy region between 3-10~ 3 and 3- 10 _1 eV. Assuming two different 
shapes for the unknown infrared photon spectrum in this energy range we calculate 
upper limits on the infrared photon density on the basis of the power daw fit obtained 
for the observed spectrum up to the maximum energy. 
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1 Introduction 



The cosmologically important extragalactical diffuse infrared background 
(DIRB) in the astronomical window from the optical to the infrared (IR) has 
not yet been directly determined experimentally due to large sytematic errors 
driven by local effects. As realized already a long time ago, the detection of ex- 
tragalactic TeV-7 sources would enable one to indirectly measure this photon 
density due to unavoidable pair production losses of TeV photons in infrared 
photon fields (Gould & Schreder 1966, Stecker et al. 1992). Interpreting the 
Whipple data of Mkn421 (redshift z = 0.031), Stecker et al. (1994) claimed 
to have seen an exponential cut-off resulting from cosmic absorption of an 
otherwise smooth power law 7-ray spectrum. The authors suggested that the 
IR-density is given by n(e) « (0.8±g;t) ■ 10-V 2 - 6 (/i/0.75) cm^eV- 1 where h 
refers to the normalization factor of the Hubble constant (H = /ilOOkm s" 1 
Mpc" 1 ). A similar analysis has been carried out by Dwek & Slavin (Dwek & 
Slavin 1994). The inferred large values of the diffuse near-infrared background 
density extrapolated to the infrared would imply that it is virtually impossi- 
ble to discover any extragalactic source above a few TeV. However, Biller et 
al. (Biller et al. 1995) correctly pointed out that unless the source spectrum 
is known, one measured TeV source can only yield an upper limit, because 
the cut-off may be due to internal absorption at the source. They obtain a 
conservative upper limit of e 2 n(e) = 0.04eVcm~ 3 at e = 0.1 eV. 
In the meantime, besides Mkn421 (Punch et al. 1992, Petry et al. 1996) 
Mkn501 (redshift z = 0.034) is now the second extragalactic TeV source 
to be discovered and extensively monitored with the Whipple and HEGRA 
Cerenkov telescopes (Kerrick et al. 1995a, Bradbury et al. 1997). Both ob- 
jects belong to the blazar subclass of galaxies showing powerful non-stellar 
activity characterized by rapidly variable, polarized continuum emission. As 
the outburst of Mkn501 in early March 1997 (Aharonian et al. 1997) and the 
two outbursts of Mkn421 in Spring of 1995 and 1996 (Kerrick et al. 1995b, 
Gaidos et al. 1996, Buckley et al. 1996) have shown, both objects also exhibit 
rapid variability with large amplitudes in the TeV energy range. The HEGRA 
collaboration showed that between March 16 and March 20 (about 27 hours 
observation time) the Mkn501 energy spectrum in the TeV energy range ex- 
tended beyond 10 TeV without a visible break in the power law spectrum with 
a differential power law index of 2.49 ± 0.11 (stat.) ± 0.25 (syst.) (Aharonian 
et al. 1997). During the quiescent state of Mkn501 in 1996, this source was 
observed with the HEGRA CT1 telescope for 220 hours. The measured energy 
spectrum for this period can be described by a differential power law index of 
2.5 ± 0.4 (total error) (Petry 1997). 

In the following we first discuss the current situation regarding models and ob- 
servations of the DIRB and then derive an upper limit on the DIRB from the 
measured energy spectrum of Mkn501 in 1997 for IR photon energies between 
3-10" 3 eV and 3-10" 1 eV. 
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2 Models and Observations of the DIRB 



The current experimental and model situation regarding the DIRB is summa- 
rized in fig. 1. Since the published energy spectra for both of the extragalactic 
sources show no spectral break we do not incorporate the claimed detections 
of the DIRB density based on preliminary Mkn421 spectra (Stecker et al. 
1994). We include, however, the recent tentative detection of the far-infrared 
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Fig. 1. Energy density of the extragalactic diffuse background radiation; stars: 
tentative FIRAS detection of the FIR background without CMBR (Puget et al. 
1996), full downward triangles: upper limits from the DIRBE experiment ([Hauser, 
p. 135], all references in square brackets are taken from Calzetti et al. 1995), squares: 
lower limit from faint blue galaxies [Tyson, p. 103], circles: upper limit of possible 
detection [Paresce, p. 307], open triangle: IRAS upper limit (Boulanger & Perault 
1988), full upward triangles: evolution model dependent lower limits from number 
counts of infrared-bright galaxies from Hacking & Soifer 1991, solid line: average 
model from MacMinn & Primack including the CMBR (MacMinn & Primack 1996), 
dashed line: average model from Fall, Chariot & Pei added to the CMBR (Fall et 
al. 1996), dotted line: upper limit derived in this paper. 

background radiation (Puget et al. 1996) by COBE which is an important 
step towards a direct measurement of the flux in the IR-regime and which is 
at the moment rather weakly constrained by the results given by Hauser in 
Calzetti et al. (1995). 
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Due to the cosmological implications of the DIRB photon density a great 
number of diffuse IR models ranging from simple power laws to multicom- 
ponent spectra have been developed by many authors. A large number of 
parameters such as the star/galaxy formation rate, the number distribution of 
star masses, the dust content and all cosmological parameters enter into these 
models. In fig. 1 we show the predictions of two complex models: The model 
which we later employ for the numerical calculation is by MacMinn & Pri- 
mack (MacMinn & Primack 1996) who have calculated the DIRB for various 
realizations of cosmological parameters and dark matter models. The results 
clearly depend on the set of parameters, and the flux predictions vary at most 
by a factor of 3. Accounting for the uncertainties of the model an average 
DIRB spectrum is assumed in the following. The model by Fall, Chariot & 
Pei (Fall et al. 1996) is based on the assumption that the star formation rate is 
directly related to the consumption of neutral gas. The density of neutral gas 
can be determined from analyses of Lya absorption lines as seen in distant 
quasar spectra thus leading to a cosmological rate of star formation out to 
distances of z ~ 5. Once the rate is known, the emissivity of the universe in 
the IR can be calculated as a function of z. The dashed curve in fig. 1 shows 
the case where the neutral gas is fully consumed during the formation process 
referred to as the closed box scenario. 



3 Gamma-ray absorption 



For 7-rays of energy E propagating from a distant source at redshift z Q towards 
a terrestrial observer the threshold energy for pair creation in interactions 
with low energy photons of present-day energy e from an isotropic diffuse 
background radiation field is given by 
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where \i = cos 9 denotes the cosine of the scattering angle. A soft photon 
density strongly varying with energy is thus reflected in the optical depth 
determining the number density of target photons at the resonant energy 
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cx E . The pair creation cross section is given by 
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Here <tt = 6.65 • 10~ 25 cm 2 denotes the Thomson cross section. For the com- 
putation of the optical depth we use the geodesic radial displacement function 
dl/dz = + z)E(z)]~ 1 . With the proper physical distance, l(t), between a 
pair of well-separated points as a function of time given by l(t) = loa(t), the 
cosmological expansion rate is H(t) = ^ = H Q E(z) with the dimensionsless 
function E(z). Under the assumption that the mean mass density is domi- 
nated by non-relativistic matter the cosmological equation for the expansion 
rate is given by (see e. g. Peebles 1993, eq. (13.3)) 

- = H E(z) = H [Q(1 + zf + Q R (1 + zf + Ov] 1/2 (4) 

(X 



where f2, Q R , and are the three contributions to the Hubble constant due 
to the present mean mass density, the radius of space curvature, and the 
cosmological constant A, respectively, with fi+f^+f^A = 1. For a cosmological 
model with Q = 1, A = 0, and negligible space curvature the function E(z) 
simplifies to (1 + z) 3 ^ 2 and the optical depth can be written as 
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for a non-evolving present-day background density rib, i-e. n' h (z, e')de' = (1 + 
z) 3 n^(e)de where the prime indicates comoving frame quantities. Numerical 
results for the optical depth using specific models for the background radia- 
tion will be obtained below. 

Under the assumption of a particular model for the IR-density one can nu- 
merically integrate equation (5) and obtain the optical depth r 77 (_E, z). For a 
source spectrum $(£") the observed spectrum is then simply given by $(£') x 
exp(— r 77 (S, z)). While for an optical depth smaller than unity the universe 
appears transparent, it becomes opaque for larger values of t^ 7 (E,z). This 
defines the 7-ray horizon (r 77 (_E, z) — 1) which is the size of the visible 7-ray 
universe at a given energy. The 7-ray horizon is shown for two values of the 
Hubble constant in fig. 2. It turns out that at energies between 10 10 and 10 15 eV 
the size of the visible universe is a decreasing function of energy. Note that 
the adopted model does not take into account the extragalactic diffuse radio 
background which only becomes important above the characteristic energies 
of air showers accessible by HEGRA-type experiments. 
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Fig. 2. 7-ray horizon for the averaged diffuse background radiation model from 
MacMinn & Primack (fig. 1) Solid line: h = 1.0. Dotted line: h = 0.5. To guide the 
eye the DIRB model has been added to the figure with an inverse energy scale. 

4 Results 

Petry (Petry 1997) and Aharonian et al. (Aharonian et al. 1997) have shown 
that the TeV energy spectrum of Mkn501 can be described by an unattenu- 
ated power law with differential indices of —2.5 ±0.4 (total error) in 1996 and 
—2.49 ± 0.11 (stat.) in May 1997, respectively. For the analysis of the 1996 
data taken with the HEGRA CT1 telescope during about 220 hours of obser- 
vation time in fig. 3 we show the integral flux spectrum with the fitted power 
law index of —1.5 ± 0.3 (total error) indicating the reduced error in the fit of 
the integral spectrum due to the correlated bins. This flux is compared to the 
Crab nebula integral spectrum as measured in the 1995/96 observation period 
with the CT1 telescope (37 hours). From the analysis of the 1997 Mkn501 (27 
hours) and Crab nebula data (10 hours) as measured with the HEGRA CT 
system we show the differential flux spectra in fig. 4. From the observed unat- 
tenuated energy spectrum extending up to 10 TeV we conclude that the optical 
depth for this source is less than unity, i.e. , r 77 (i? = 10 TeV, z = 0.034) < 1. 
To derive the upper limit on the DIRB photon density we take r 77 = 1.0 and 
h = 1.0, and assume a specific shape of the IR-density. 

For the energy range 3 • 10 _3 eV < e < 3 • 10 _1 eV we make an ansatz for the 
infrared spectrum based on the results of two different models (outside this 
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Fig. 3. Average integral spectrum of 7-rays from Mkn 501 as measured in 1996 (220 
hours) and from the Crab nebula as measured in the period 1995/96 (37 hours) 
with the HEGRA CT1 telescope. The lines represent power-law fits. The error bars 
correspond to the statistical errors (from Petry 1997). 




Fig. 4. Average differential spectrum of 7-rays from Mkn 501 between March 15 
and March 20, 1997 (27 hours) and from the Crab nebula (10 hours) as measured 
with the HEGRA CT system. The Crab data points are scaled by a factor of 0.2. 
The lines represent power-law fits. Only statistical errors are shown. The energy 
scale has a 20% systematic error (from Aharonian et al. 1997). 
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Fig. 5. Results from this analysis: Values for the spectral index v and n above the 
upper curve are forbidden given the unattenuated observation of Mkn 501 at 10 TeV 
(squares). For illustration purposes we add the exclusion curve that would result if 
no spectral break would be observed up to 15 TeV (circles). 

energy range only the MacMinn & Primack model shown in fig. 1 is used): (i) 
e 2 n(e) = n (t/E p ) v with the pivot energy point E p = 3 • lCT 2 eV and (ii) the 
shape of the spectrum in this energy range is assumed to be identical with the 
models by MacMinn & Primack but its absolute level is varied in the calcula- 
tion. The results for the first ansatz are presented in fig. 5 where values of v 
and n Q below the curve are still allowed by the detection of Mkn 501 at 10 TeV. 
To illustrate the impact of a still not excluded observation of the unabsorbed 
spectrum extending up to 15 TeV, a second curve was added in the figure. 
Using the second model ansatz the observation of Mkn 501 at 10 TeV restricts 
the normalization N of the spectral segment between 3 • 10 -3 eV and 3 • 10 _1 eV 
to N < 1.8 (times the flux assumed in this work). At the pivot energy point 
of 3 • 10~ 2 eV we thus determine upper limits of e 2 n(e) = 1.1 • 10~ 3 eV/cm 3 for 
the MacMinn & Primack model and 1.0 • 10~ 3 eV/ cm 3 for the powerlaw ansatz 
and assuming the slope of the spectrum to be flat (v = 0) around the pivot 
point. A comparison with fig. 1, where we have added the upper limit curve 
derived using the MacMinn & Primack ansatz, indicates that these limits are 
about two orders of magnitude below those accessible by direct measurements 
at this energy and that they are compatible with the tentative FIRAS mea- 
surement of the IR density (Puget et al. 1996). 

As pointed out in the introduction, the measurement of single spectra, even 
if a break-off feature is observed, does not permit a determination of the 
actual level of the DIRB density. Up to the end of the 1997 observation pe- 
riod (October 1997) the blazar Mkn 501 continued to be in a high state with 
respect to TeV emission. The increasing statistics will enable Cerenkov tele- 
scope experiments to extend observations to higher energies and search for 
cut-off features. The relationship between the infrared photon density in the 
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energy range from 3 • 10~ 3 eV to 3 • 10 _1 eV, assuming a powerlaw spectrum 
e 2 n(e) = n (e/3 • 10~ 2 eV)~ a4 , and the maximum observable energy, i. e., E max 
with optical depth r = 1, for different values of redshift z is shown in fig. 6. 
If e. g. no cut-off is observed up to around 70 TeV we would thus derive, using 
the above described procedure, an upper limit on the DIRB density of about 
0.1 times the average MacMinn & Primack model prediction. 
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Fig. 6. The normalization uq of the infrared photon density in the energy range 
between 3 • 1CT 3 to 3 • 1CT 1 eV assuming a spectrum e 2 n(e) = no(e/3 • 1CT 2 eV)~ a4 
which would result from our calculations as a function of the maximum observable 
energy (with optical depth r = 1) for three values of redshift z. 



5 Conclusions 

Based on the observation of the unabsorbed Mkn501 y-spectrum extending 
beyond 10 TeV, we have derived a stringent upper limit on the extragalactic 
diffuse infrared photon energy density in the energy range from 3 • 10~ 3 to 
3 • 10 _1 eV of 1.8 times the prediction of an average model by MacMinn & 
Primack. This translates into an upper limit of the energy density at the used 
pivot energy point of 3 • 10~ 2 eV of e 2 n(e) = 1.1 • 10~ 3 eV/cm 3 . For the sec- 
ond ansatz, i. e. , a power-law ansatz for the DIRB around this pivot energy 
point, we determined upper limits on the normalization as a function of the 
spectral index, e.g. for a flat spectrum e 2 n(e) = Hq, the resulting upper limit 
is n = 1.0 • 10 _3 eV/cm 3 . These upper limits are about 2 orders of magnitude 
below upper limits derived from current direct measurements at this energy 
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and are not in contradiction to a value of the DIRB density derived from pre- 
liminary evidence of TeV emission reported by Meyer & Westerhoff ( Meyer 
& Westerhoff 1996). 

The results presented in this paper are well compatible with other analyses 
of the infrared photon density based on the HEGRA CT data of Mkn501 
(Mannheim 1997, Stanev & Franceschini 1997). They are also in good agree- 
ment with an empirical calculation of the DIRB based on galaxy luminosity 
functions in the IR (Malkan & Stecker 1997) which in turn is in good agree- 
ment with recent DIRB models derived from star formation data (Guiderdoni 
et al. 1997). 
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